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Nonlocality and quantum entanglement constitute two special aspects of the quantum correlations 
existing in quantum systems, which are of paramount importance in quantum-information theory. 
Traditionally, they have been regarded as identical (equivalent, in fact, for pure two qubit states, 
that is, Gisin's Theorem), yet they constitute different resources. Describing nonlocality by means 
of the violation of several Bell inequalities, we obtain by direct optimization those states of two 
qubits that maximally violate a Bell inequality, in terms of their degree of mixture as measured by 
either their participation ratio R = 1/Tr(p 2 ) or their maximum eigenvalue \ ma x- This optimum 
value is obtained as well, which coincides with previous results. Comparison with entanglement is 
performed too. An example of an application is given in the XY model. In this novel approximation, 
we also concentrate on the nonlocality for linear combinations of pure states of two qubits, providing 
a closed form for their maximal nonlocality measure. The case of Bell diagonal mixed states of two 
qubits is also extensively studied. Special attention concerning the connection between nonlocality 
and entanglement for mixed states of two qubits is paid to the so called maximally entangled mixed 
states. Additional aspects for the case of two qubits are also described in detail. Since we deal 
with qubit systems, we will perform an analogous study for three qubits, employing similar tools. 
Relation between distillability and nonlocality is explored quantitatively for the whole space of 
states of three qubits. We finally extend our analysis to four qubit systems, where nonlocality for 
generalized Greenberger-Horne-Zeilinger states of arbitrary number of parties is computed. 

PACS numbers: 03.65.Ud; 03.67.Bg; 03.67.Mn 



I. INTRODUCTION 

Schrodinger's [1] modern notion of entangled state 
historically appeared within the debate of the paradox 
posed by Einstein, Podolsky and Rosen (EPR) [2]. EPR 
pointed out the possible lack of completeness of the new- 
born theory of quantum mechanics. In their famous pa- 
per they suggested a description of nature (called "local 
realism") which assigns an independent and objective re- 
ality to the properties of separated parties of a composite 
physical system. EPR applied the criterion of local real- 
ism to predictions associated with an entangled state, a 
state that cannot be described solely in terms of the prop- 
erties of its subsystems, to conclude that quantum me- 
chanics was incomplete. Schrddinger, instead, regarded 
entanglement as the characteristic feature of quantum 
mechanics. Quantitatively though, no measure for this 
"quantum strangeness" was provided at the time. 

The most significant progress toward the resolution of 
the EPR debate came with Bell's work. Bell [3, 4] proved 
the impossibility of reproducing all correlations observed 
in composite quantum systems using models similar to 
that of EPR. In fact, Bell showed that local realism, in 
the form of local variable models (LVM), implies con- 
straints on the predictions of spin correlations in the form 
of inequalities, also known as Bell's inequalities, which 
can be violated by quantum mechanics. That is why 
quantum mechanics is regarded as being inherently non- 
local. 
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If quantum mechanics could be described by LVM, 
then the correlation values measured between parties 
could be reproduced assuming the corresponding op- 
erators had already a definite value previous to mea- 
surement. Let us consider the case of the two-party 
two-outcome scenario. Two distant observers possessing 
a distributed state pab exploit the correlations arising 
from 

P(a, b\A, B; pab) = Tr[p AB (U^ ® nf )] (1) 

quantum mechanically, and are positive oper- 
ators involving measurements of observables A and B, 
with outcomes a and b satisfying AU^ = all^ and 
£?nf = fellf . If a LVM could mimic the same corre- 
lations 

P(a, b\A, B; p AB ) = J2 P W A < X ) P ( b \ B > A )-"(A), (2) 

A 

with n(X) being a probability measure for the classical 
variable A and P(a\A,X) and P(b\B, A) local functions, 
therefore the use of state pab would provide no improve- 
ment over classical resources. Thus it was clear that the 
notion of nonlocality of a state pab would emerge if there 
existed no LVM that through (2) could reproduce the 
quantum mechanical results of (1). Nonlocality was a 
character of entangled states via violation of a Bell in- 
equality. 

Ever since Bell's contribution, entanglement and non- 
locality became similar terms. The nonlocal character of 
entangled states was clear for pure states. In fact, all 
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entangled pure states of two qubits violate the CHSH 
inequality and therefore are nonlocal [5]. The situation 
became more involved when Werner [6] discovered that 
while entanglement is necessary for a state to be nonlo- 
cal, for mixed states is not sufficient. He also introduced 
the usual (modern) definition of entangled state: given 
two parties A and B, a shared state pab is termed un- 
entangled or separable if it cannot be expressed as the 
mixture of product states 

JV 

pab = J2 ^Iv4><v^l®iv4><v4l, (3) 

i=l 

that is, when its preparation does not require a nonlocal 
quantum resource. This definition, in spite of its clear 
physical meaning, is somewhat impractical, since tests 
to distinguish separable from entangled states are com- 
plicated. 

With the advent of quantum-information theory 
(QIT), the interest in entanglement has dramatically in- 
creased over the years since it lies at the basis of some 
of the most important processes and applications studied 
by QIT such as quantum cryptographic key distribution 
[7], quantum telcportation [8], supcrdense coding [9] and 
quantum computation [10, 11], among many others which 
possess no classical counterpart. All these tasks require 
distributed quantum correlations between parties, and 
the only means available in nature are entangled states. 
Using a modern nomenclature, when a quantum state 
cannot be prepared using only local operations and clas- 
sical communications (LOCC), it is said that it possesses 
quantum correlations and the state is entangled. Spa- 
tially separated observers sharing an entangled state and 
performing measurements on them may induce nonlocal 
correlations which cannot be simulated by local means 
(violate Bell inequalities). 

Confusion between nonlocality and entanglement ap- 
peared during the finer study of the usefulness of quan- 
tum correlations. Entanglement is commonly viewed as a 
useful resource for various information-processing tasks. 
Yet, there exist certain tasks, such as device-independent 
quantum key distribution [12] and quantum communica- 
tion complexity problems [13], which can only be carried 
out provided the corresponding entangled states exhibit 
nonlocal correlations. Then we are naturally led to the 
question of whether nonlocality and entanglement con- 
stitute two different resources. 

The purpose of the present work is to shed some light 
upon the relation between entanglement and nonlocality 
through the maximal violation of a Bell inequality for 
two, three and four qubit systems. Throughout the ar- 
ticle, and in order to avoid confusion, we will refer to 
the quantity "nonlocality" as being equivalent to "maxi- 
mum violation of a Bell inequality" . However, the usual 
meaning of nonlocality (or that of a nonlocal state) as a 
concept involving the mere violation of a Bell inequality 
remains the same. Although detection and characteriza- 
tion of entanglement is far from being complete, its status 



is more developed than that of nonlocality. The problems 
experienced in defining a unique measure of entanglement 
in the multipartite case, in the form of partitions in the 
system, disappear in the nonlocality case since there ex- 
ists well defined Bell inequalities for multipartite systems 
of arbitrary number of qubits. This fact makes the study 
of nonlocality conceptually and quantitatively easier. 

This paper is organized as follows. In Section II we 
review recent results concerning nonlocality in bipartite 
physical systems and concentrate on the CHSH Bell in- 
equality [14] for two qubits. We also obtain, after a di- 
rect optimization over the observers' settings, the family 
of two qubit mixed states that optimizes the violation of 
the CHSH inequality for a given degree of mixture, as well 
as the concomitant optimal value for CHSH. This novel 
approach recovers and extends previous results found in 
the literature. Our results concerning the connection be- 
tween nonlocality and entanglement for mixed states of 
two qubits also pay special attention to the case of the 
so called maximally entangled mixed states (MEMS). As 
far as the duality entanglement-nonlocality is concerned, 
these very same results find too interesting echoes in a 
well known condensed matter system, namely, the infi- 
nite XY model. How nonlocality can be present in lin- 
ear combinations of pure states of two qubits is also re- 
ported, where the optimal value of the violation of the 
CHSH inequality is obtained for any given superposition. 
Parallelism with the role of entanglement of superposi- 
tion is also discussed. Section III is devoted to the study 
of entanglement and nonlocality for three qubit systems, 
employing similar tools. Relation between distillability 
and nonlocality is explored quantitatively for the whole 
space of states of three qubits. Section IV extends the 
present subject of study to four qubit systems, where 
nonlocality for generalized Greenberger-Horne-Zeilingcr 
(GHZ) states of arbitrary number of parties is computed. 
Finally, some conclusions are drawn in Section V. 



II. TWO QUBITS 
A. Nonlocality and the CHSH Bell inequality 

LVM cannot exhibit arbitrary correlations. Mathemat- 
ically, the conditions these correlations must obey can 
always be written as inequalities -the Bell inequalities- 
satisfied for the joint probabilities of outcomes. We say 
that a quantum state p is nonlocal if and only if there 
are measurements on p that produce a correlation that 
violates a Bell inequality. 

Most of our knowledge on Bell inequalities and their 
quantum mechanical violation is based on the CHSH 
inequality [14]. With two dichotomic observables per 
party, it is the simplest [15] (up to local symmetries) 
nontrivial Bell inequality for the bipartite case with bi- 
nary inputs and outcomes. Let A\ and A 2 be two 
possible measurements on A side whose outcomes are 
aj e {— 1,+1}, and similarly for the B side. Math- 
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ematically, it can be shown that, following LVM (2), 
\^chshW\ = Wih + «i & 2 + a 2 h - a 2 b 2 \ < 2. Since 
ai(&i) and a 2 (b 2 ) cannot be measured simultaneously, 
instead one estimates after randomly chosen measure- 
ments the average value Bggg H = B ohsh( x )pW = 
E{A U B X ) + E(A U B 2 ) + E(A 2 ,B 1 ) - E(A 2 ,B 2 ), where 
E(-) represents the expectation value. Therefore the 
CHSH inequality reduces to 



Kiwi < 2- (4) 

Quantum mechanically, since we are dealing with 
qubits, these observables reduce to Aj(Bj) = aj(bj) • a, 
where aj(bj) are unit vectors in R 3 and a — (a x ,a y , <r z ) 
are the usual Pauli matrices. Therefore the quantal pre- 
diction for (4) reduces to the expectation value of the 
operator B C hsh 



B. Maximal violation of the CHSH inequality and 
mixedness for two qubits 

What is the maximum violation of Bchsh{= 
Tr(pBcHSH) < 2) for a given state p? 

Before any attempt to proceed with a definite optimiza- 
tion program, we shall undertake a detailed analysis of 
the special form for Bqhsh — Tr(pBcHSH), the basic 
quantity we are about to deal with. The nature of any 
bipartite mixed state of two qubits is described by a pos- 
itive, semi definite matrix, whose eigenvalues {A^} are 
such that they < Xi < 1 and ^ A, = 1. In other 
words, the complete description of the density matrix p 
necessitates 4 2 — 1 = 15 real parameters. 

Usually, the preferred basis for two qubit states is the 
so called computational basis {|00), |01), |10), |11)}. In 
our case, it will prove convenient to employ the so called 
Bell basis of maximally correlated states, which are of 
the form 



Bi 
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Bi - A 2 



B 2 



(5) 



Tsirelson showed [16] that CHSH inequality (4) is maxi- 
mally violated by a multiplicative factor \[2 (Tsirelson's 
bound) on the basis of quantum mechanics. In fact, it 
is true that \Tr(pABBcHSH)\ < 2-\/2 for all observables 
Ai, A 2 , Bi, B 2 , and all states pab- Increasing the size 
of Hilbert spaces on either A and B sides would not give 
any advantage in the violation of the CHSH inequalities. 
In general, it is not known how to calculate the best such 
bound for an arbitrary Bell inequality, although several 
techniques have been developed [17]. 

A good witness of useful correlations is, in many cases, 
the violation of a Bell inequality by a quantum state. 
But not all entangled states are nonlocal. Although this 
is the case for pure states of two qubits (CHSH inequality 
violation), Werner showed that it cannot be generalized 
to mixed states. After introducing the states which are 
now called Werner states 
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(7) 



Now we rise the question of wether all elements of p in- 
tervene in the computation of Bqhsh- Given a state p, 
we make a change of basis so that we work in the Bell 
basis. For simplicity, and without loss of generality, we 
shall take real coefficients (9 = in (7)). Due to the 
special form for (5), we separate the elements of p (now 
in the Bell basis) into two contributions, namely 



P = P\\ + P± 
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(9) 



pw =f# I + (1 -P) : 



(6) 



where \il>~) is the singlet state and / is the 4x4 identity, 
he provided a LVM for measurement outcomes for some 
entangled states of this family. Although promising new 
results have been obtained recently [18, 19], even in the 
simplest case of Werner states of two qubits (6), it is 
in general extremely difficult to determine whether an 
entangled state has a LVM or not, since finding all Bell 
inequalities is a computationally hard problem [20, 21]. 

Therefore we shall consider the optimization of the vio- 
lation of the CHSH inequality over the observer's settings 
as a definitive measure for both signaling and quantifying 
nonlocality in two qubit systems. 



This separation is motivated by the fact that only terms 
in p|| contribute to Tr(pBcHSH), as one can easily 
check. In other words, Tr(pBcHSH) — Tt{p\\Bchsh) + 
Tr(p ± B CHSH ) = Tr{p\\B CHSH ). The superscripts of 
the matrix elements in (8) refer to the concomitant real 
(R) and imaginary (/) parts. This observation consti- 
tutes the starting point of our study. 

The answer to the initial question of this Section in- 
volves only the elements of py for p in (8). The latter fact 
enormously simplifies the computation, but we neverthe- 
less encounter a highly nontrivial optimization enterprise. 

Fortunately, we do not require all elements of p\\ . In- 
stead, since we seek maximum nonlocality, we will con- 
sider states which are diagonal in the Bell basis (null 
elements off-diagonal in p\\ (8)), for nonlocal correlations 
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concentrate after some depolarizing process [22]. Previ- 
ous authors computed the entanglement and the maxi- 
mum violation for the Bell inequality for Bell diagonal 
states [23], and the particular form for these states [24]. 
In the present work, we re obtain and extend their results 
by means of a specific optimization technique, which is 
described in full detail in Appendix I. 
For diagonal states in the Bell basis 



the latter result with relation (14), we finally [25] arrive 
at 

B^hsh = ^Tr[Bl HSH ] ■ = 4 ■ yj^£. (15) 

This result is valid for the range R e [2,4]. The corre- 
sponding state (13) can now be cast in the new form 



£ S) =Ai|$ + )($ + | + A 2 |$-)($-| 

+ A 3 |*+)(* + | + A 4 |*-)(*-|, (10) 

with eigenvalues appearing in decreasing order, we obtain 

m a xTr(pf e fBcHSH) = 2y/2y/(X 1 - A 4 ) 2 + (A 2 - A 3 ) 2 . 

(11) 

Recall that 2y/2 is the maximum value allowed by quan- 
tum mechanics (attained only for states (7)). 

We are going to determine which is the maximum ex- 
pectation value of the CHSH operator (5) that a two 
qubit mixed state p with some degree of mixedness, 
in this case given by the so called participation ratio 
R = 1/Tr(p 2 ), may have. Notice that no assumption 
is needed regarding the state being diagonal or not in 
the Bell basis. In order to solve the concomitant vari- 
ational problem (and bearing in mind that Bchsh — 
Tr(pBcHSH)), let us first find the state that extrem- 
izes Tr(p 2 ) under the constraints associated with a given 
value of Bchsh, and the normalization of p. This vari- 
ational problem can be cast as 



S[Tr(p 2 ) + (3Tr(pB C HSH) - aTr(p)] = 0, (12) 

where a and /3 are appropriate Lagrange multipliers. 
The solution of the above variational equation is given 

by 



P = 2 [ aI ~ P®chsh] ■ 



(13) 



with / being the 4x4 identity matrix. The value of the 
Lagrange multiplier a is immediately obtained by the 
normalization requirement, with a = |. 

From (13) we find Bchsh, multiply it by p and ap- 
ply the corresponding definition of £?™ff SH , taking into 
account that Tr(p 2 ) = 1/R. We arrive at the result 



Bchsh = Tr(pB C HS H ) = - | ■ (14) 



By either squaring (13) and taking the trace according 
to the definition of R, or rather multiply it by Bchsh 
(Bchsh is traceless) in order to get B C hsh, both ways 
lead to f3 = — Bc " SH , Combining either the former or 



p u = diag^x,x, 1 — ^ 



2x l-2x 



(16) 



with x e [0, |] and diagonal in the Bell basis. In the 
region R e [1,2] the form of state (16) is no longer valid. 
Instead, we look for those states that stay close to pure 
states (maximum nonlocality) and possess rank 2 (fol- 
lowing the requirement Re [1,2]). The simplest case is 
that of a state diagonal in the Bell basis, being of the 
form 



Pi = diag(l — x, x, 0, 0) 
with x e [0, h]. Relation (11) returns 



B^sh(pi) - 2V2^(l-x) 2 + x 2 



(17) 



(18) 



We are now in a position to answer the initial question 
of this Section. Besides, we do not find the functional 
form for B^ff SH (R) for two qubits, we do also obtain 
the form for those states, states which are diagonal in the 
Bell basis. We shall call these states Maximally Nonlocal 
Mixed States (MNMS). 

For a given value of the participation ratio R, we can 
obtain a more general class of states rather than pj (17) 
and pu (16) by letting a non-zero phase 9 in (7). By do- 
ing so, and rewriting the concomitant MNMS states pj 
(17) in the computational basis, one easily obtains the 
family of states provided in [24]. But not only this: a 
whole series of new states possessing maximum nonlocal- 
ity for a given value of R are obtained by changing the 
position of the eigenstates of pi in (17). This is possible 
since no preferred disposition of states in the Bell basis 
is required for diagonal states (35) as far as maximum 
amount of nonlocality is concerned. 

Let us summarize all previous results: the maximum 
amount of nonlocality attained by a mixed state p of two 
qubits is given by 

• Mixedness described by the participation ratio R = 
1/Tr(p 2 ) 



°CHSH\ n ) 



B^ SH {R)=A- 



Me [1,2] 



I A- R 
~4R~ : 



(19) 
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• Mixedness described by the maximum eigenvalue 

Amax(p) 
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(20) 



One must bear in mind that states that reach the 
maximum possible value for nonlocality measure B^^ SH 
greatly depend on what measure for the degree of mix- 
ture is employed. The only case where both descriptions 
agree is for those states that strictly violate the CHSH 
inequality, namely, the MNMS pi (17). 



C. Nonlocality for maximally entangled mixed 
states 

Maximally entangled mixed states (MEMS) constitute 
a family of states that are maximally entangled for a 
given degree of mixture, measured by the participation 
ratio R = 1/Tr(p 2 ). In practice, one will more often 
have to deal with mixed states than with pure ones. 
From the point of view of entanglement-exploitation, one 
should then be interested in MEMS states pmems, which 
are basic constituents of quantum communication proto- 
cols. The MEMS states have been studied, for example, 
in Refs. [26-28]. MEMS states have been experimen- 
tally encountered [29, 30]. In the computational basis 
{|00), |01), 1 10) , |11)}, they are written as 




x/2 

2g(x) 



g{x) 



(21) 



with g(x) = 1/3 for < x < 2/3, and g(x) = x/2 for 
2/3 < x < 1. The quantity x is equal to the concurrence 
C. The change of g(x)— regime ensues for R = 1.8. 

Our goal is to uncover interesting correlations between 
entanglement, nonlocality and mixedness that emerge 
for these states. Indeed, the study of the nonlocality 
of these states offers an excellent framework where to 
compare the extremal cases for nonlocality and entan- 
glement. The MEMS states (21) are written in the Bell 
basis {|$ + ), |<&~), |* + ), |* - )} in the form 



/ g(x) ^ 
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(22) 



which is to be compared with the general form for arbi- 
trary states (8) in the Bell basis. The direct comparison 
of MEMS states in the Bell basis yields to the conclusion 
that states Pmems (22) behave as if they were diagonal 
in the Bell basis as far as nonlocality is concerned. 

This crucial observation allow us to the simple calcu- 
lation of the maximum amount of nonlocality for MEMS 
states to be of the form 



omax _ ffvT+9i* 0<*<i 

Bchsh[x) ~ \2V2x, §<*<i 



(23) 



Recall that x = C, the so called concurrence. Also, it is 
plain from relation (23) that any bipartite state possess- 
ing x > will violate the CHSH inequality, since no 
state is more entangled that the MEMS states. 

The concurrence entanglement indicator C for MNMS 
states pi (17) can be easily calculated to be C = 1 — 
2x. Therefore, the relation between nonlocality B™ff SH 
and C is such that B™ff SH = 2\f\ + C 2 , which easily 
recovers a previous result [23]. In other words, for a given 
value of C, MNMS states possess the maximum possible 
violation of the CHSH inequality while, on the contrary, 
MEMS possess a minimum amount of CHSH violation 
(2V1 + C 2 > 2V2C). 

As a consequence, we draw the conclusion that max- 
imum entanglement for mixed states of two qubits does 
not imply maximum nonlocality, though in this extremal 
case nonlocality and entanglement are monotonic in- 
creasing functions one of the other. 



D. A physically-motivated case: the XY model 

The general study of nonlocality and entanglement in 
an infinite quantum system was performed in Ref. [31]. 
In this section we incorporate new results and generelize 
previous ones in the light of the bounds encountered for 
the maximum violation of the CHSH Bell inequality for 
a given degree of mixture. 

The general two-site density matrix for two spins along 
the XY chain is expressed as 



p {R) = 1 



/ j UV ^ 11 



(24) 



R = j — i indicates the distance between spins (not to 
be confused with the participation ratio R = 1/Tr(p 2 ), 
{u,v} denote any index of {er , a x , cr y , c z }, and Tuv = 
{a l u ® a 3 v ). Due to symmetry considerations, only 
{T^\T^y\T^\T^y ] } do not vanish. Barouch et al [32] 
provided exact expressions for two-point correlations, to- 
gether with all the dynamics associated with an external 
magnetic field h(t) along the z-axis. We shall consider 
the case where h jumps from and initial value ho to a 
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final value hf at t — 0, that is, a quench (the equilibrium 
case is easily recovered when hf = ho) and the R = 1 
configuration (nearest neighbors). 

The most remarkable result of Rcf. [31] as far as non- 
locality is concerned is that the maximum value for the 
quantity B™ff SH for states (24), given by twice the ex- 
pression 

V /||TW|P-min([rW] 2 , [T™]\ [rj*>] 2 ) + 2 [T™\ \ 

(25) 

with T( R ) = (Tj^T^.Tj?)), is always < 2 for any 
configuration R and any non-zero value of the entangle- 
ment, 

Fig. 1 depicts several time evolution for the state 
(24) once a quench in the external magnetic field is ap- 
plied. As a consequence, we have a nonergodic evolution 
in time, which translates into oscillating values for both 
Bqjj SH and R. The previous time dependent cases cor- 
respond to Fig. 1(a) and Fig. 1(b). The static case is de- 
picted in Fig. 1(c), where nonlocality improves for states 
approaching the Ising case (7 = 1). Several time evolu- 
tion plots appear in Fig. 1(d) for B™^ SH and the entan- 
glement of formation for states (24), together with the 
thermodynamic magnetization M z after a quench from 
ho = 0.5 to hf — 0. All these three quantities possess 
a nonergodic behavior, which is not surprising for they 
ultimately depend on two point spin correlators, which 
in turn are nonergodic quantities [31]. 

In Fig. 2, we consider the maximum value for Bq^j SH 
and the concurrence C that any state of the type (24) can 
have. The value for the factorizing field h s = yl — 7 2 
for which states (24) are separable [32] correspond to 
zero concurrence, that is, a line at the bottom of the 
plot. In the language of mixedness, the magnetic field 
h and the mixture of the state go in opposite directions: 
the greater the former, the lesser the latter. This is so 
because as h — > 00, we approach a pure state (R = 1) 
with both spins down. In the surface 7 = constant we 
encounter that both B^^ SH and C diverge (their first 
derivative with respect to h), thus signaling a quantum 
phase transition (except for the isotropic case 7 = 0). 

As we can see, the notion of nonlocality and entan- 
glement as different resources appear in a real physical 
system. The necessity of the violation of some Bell in- 
equality for some information-theoretic tasks and its rela- 
tion to entanglement makes the whole picture a bit more 
intriguing with this physical case. 



E. Nonlocality vs entanglement for linear 
combinations of pure states of two qubits 

Nonlocality may also exhibit interesting features for 
pure states. In our case we shall consider states of the 
form 
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FIG. 1. (Color online) (a) Value of B^ SH vs R for sev- 
eral time evolutions of the external perpendicular magnetic 
field in the XY model. The solid curve corresponds to 
the case (ho = 0.5, hf — 0). From left to right, both 
BctFsH an d R oscillate around their concomitant final (non- 
equilibrium) values. The lower dashed curve depicts the case 
(ho — 0.75, hf — 0), with a similar behavior. Notice that in 
either cases the maximum value of Bq^ sh (R) (upper solid 
line) is never crossed, (b) Similar curves for (ho = l,hf = 0) 
(solid curve) and (ho = 2, hf = 0) (lower dashed curve), (c) 
BctfsH vs R plots for several anisotropy values (from bot- 
tom to top) 7 = 0,0.1,0.3,0.5,1. As the magnetic field h 
increases from to oo (no time evolution), the curves go from 
right to left (decreasing R- values for the two qubit states). 
It is plain from this series of plots that no violation of the 
CHSH Bell inequality occurs, (d) Time evolution of BcWsh 
(lower solid curve), entanglement E (long dashed curve) and 
the magnetization M z (short dashed curve) after the quench 
(ho = 0.5, hf =0). E and M z have been shifted two units 
upwards. All these three quantities are nonergodic. See text 
for details. 



\<t>) = v / Ai|* + )+v / A^I^)+V / A^I* + )+V / Al|^>, (26) 

with real coefficients {\/A7}. This particular form clearly 
constitutes an extension of the analysis performed for 
two qubit mixed states diagonal in the Bell basis. All 
the details of the optimization of the CHSH inequality 
for states (26) is given in Appendix I. The final result is 

2V2 v /(Ai + A 4 ) 2 + (A 2 + A 3 ) 2 . (27) 

We know by virtue of Gisin's theorem [5] that pure 
bipartite entanglement implies violation of the CHSH 
inequality. The question is wether that dependency 
changes when we consider the case of linear combina- 
tions. The corresponding answer is no, because otherwise 
it would imply a preferred choice for the state basis. 

The value (27) for maximum violation of the CHSH 
inequality is to be compared with the measure of entan- 
glement for the state (26) 
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FIG. 2. (Color online) Plot of the nonlocality measure B™I?SH 
(upper surface) and twice the concurrence C (lower surface) 
for any value of the anisotropy 7 and the external perpendic- 
ular magnetic field h for two qubit states (nearest neighbors) 
in the XY model. It is plain that no violation of the CHSH 
Bell inequality occurs, regardless of the non-zero value for the 
entanglement indicator C. See text for details. 



versa) for pure two qubit states, but nothing is said re- 
garding their particular characterization. Let us illus- 
trate, before embarking on our study, what happens when 
we consider the nonlocality present in the superposed 
state 



|0>=a|Ol)+0|*+>, (30) 

with SgfWIOl}) = 2 and Bgi? SJr (|#+» = 2^2. Pre- 
sumedly, nonlocality of state \9) should be lowered by the 
action of non correlated |01). Indeed, we have 

BchshW)) = 2V2-a 2 (2-a 2 ) = 2V2-V2a 2 +0(a 4 ). 

(31) 

This example shows that BQff SH and the fidelity between 
states within the linear combination simultaneously and 
continuously change, a fact that does not occur for the 
entanglement of linear combination of states [35] . 

Theorem 1: Given a set of M orthogonal pure states 
of two qubits {\<t>i)}$L x , with YhLi <*? = 1 (<*i € R), the 
concomitant maximum nonlocality measure obeys 



C 2 = 4 det PA/B = 1 - 4(Ai + A 4 )(A 2 + A3). (28) 

Strictly speaking, the concurrence or its squared value 
C 2 are not proper measures of entanglement, for they do 
not comply with the usual basic requirements [33] . How- 
ever, they are widely used as useful entanglement quan- 
tifies (they are monotonic functions of the entanglement 
of formation Ef(-) [34], which is a good measure). 

Taking into account the form for C 2 (28), and the value 
for B™ff SH (27), we derive the relation 

B^hsh = V4C 2 +4, (29) 

which is the same result we would obtain for a pure state 
of two qubits written in the Schmidt basis. Had we as- 
sumed result (29) to hold for any state, we would have 
obtained the relation (27) for the maximum violation of 
the CHSH inequality without recourse to any optimiza- 
tion technique. 

Formula (27) has interesting echoes when compared 
to the entanglement of the same superposition of states. 
Entanglement of superposition of states was originally 
conceived in Rcf. [35], where interesting bounds for the 
superposed state were obtained in terms of their con- 
stituents. In our case, result (27) permits us to establish 
similar bounds for maximum violation of the CHSH in- 
equality. This framework offers a link between the char- 
acterization of nonlocality and entanglement, where their 
mutual intricacies become more apparent. 

As obtained before, we know by virtue of Gisin's the- 
orem [5] that nonlocality implies entanglement (and vice 



, M v M 

bswsh E «*i&> ) ^ J2( a ^ B cHSH 2 m)- (32) 

Proof: Spanning the set of states {\4>i)}f£i in the Bell 
basis, by recourse to (27) and expanding quadratic terms, 
the remaining part on the right hand side of is a strictly 
positive quantity, from hence we directly compute the 
nonlocality of superposition as stated in Theorem 1. The 

upper bound for B™ff SH 2 ( J2iLi a i|&)^ i s easily ob- 
tained by individually optimizing each term in the argu- 
ment {|</>i)}, and taking into account each contribution 
arising from 2Re[B>§% SH {{<f> i \f j ))\ , Vi f J. 

Theorem 1 and the concomitant upper bound connect 
the way nonlocality of a superposed state is distributed 
among its constituents, in a similar fashion as entangle- 
ment in Ref. [35] . More details on superposition of states 
will be described elsewhere [36]. 



III. THREE QUBITS 

A. Nonlocality for three qubit states. Application 
to the XY model 

We shall explore nonlocality in the three qubit case 
through the violation of the Mermin inequality [37] . This 
inequality was conceived originally in order to detect 
genuine three-party quantum correlations impossible to 
reproduce via LVMs. The Mermin inequality reads as 
Tr(pB M ermin) < 2, where Buermin is the Mermin oper- 
ator 
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B 



Mermin 



B, 



a 1 a 2 a 3 



B ai 



6263 



B bl 



12^3 



B hl 



b 2 a 3 1 



(33) 



with B uvw = u-cr®v-cr®w-CT with a — (cr x , <j y , a z ) being 
the usual Pauli matrices, and aj and bj unit vectors in M 3 . 
Notice that the Mermin inequality is maximally violated 
by Greenberger-Horne-Zeilinger (GHZ) states. As in the 
bipartite case, we shall define the following quantity 



We can encounter too interesting nonlocality features 
if we focus our attention to the case of pure states of 
three qubits being linear combinations in the Mermin 
basis I**) = ± 1 7 - j))/\/2- That is, states of the 
form 



(37) 



i^) = EA ± i* J ± ) 



Mermin ma * = max Tr(pB Me rmin) (34) 

as a measure for the nonlocality of the state p. While 
in the bipartite the CHSH inequality was the strongest 
possible one, this is not the case for three qubits. The 
Mermin inequality is not the only existing Bell inequality 
for three qubits, but it constitutes a simple generalization 
of the CHSH one to the tripartite case. Therefore, it will 
suffice to use this particular inequality to illustrate the 
basic results of the present work. 

In view of the previous definitions, we are naturally led 
to the question of what class of three qubit mixed states 
possesses a maximum amount of nonlocality (34), how 
does it look like and what it amounts for. Let us recall 
that the family of pure states of three qubits I**) = 

± |7 — j))/V2 forms a basis, the so called GHZ basis 
or Mermin basis, and that these states maximally violate 
the Mermin inequality. But what is the state of affairs for 
the general, mixed case? Given a state p, it can always 
be transformed into the state 



(diag) 
r Mermin 



a state which is diagonal in the Mermin basis, for non- 
local correlations concentrate after the action of some 
depolarizing process [22]. Without loss of generality, we 
can assume the eigenvalues of (35) to be sorted in de- 



creasing order, that is, Aq > A > • • • > A3 > A 3 , 



since otherwise it could be adjusted by a local unitary 
operation. 

The details of the optimization are given in Appendix 
II. However, the maximum violation of the Mermin in- 
equality is given by the quantity 



Mermin 



<4 



\ 



(36) 



with real coefficients {y Aj } such that J2j=i = 1- 

Detection and characterization of entanglement in mul- 
tipartite systems constitutes a hot research topic in QIT. 
However, no necessary and sufficient criterion is avail- 
able to date that discriminates wether a given state of a 
multipartite system is entangled or not. Indeed, highly 
entangled multipartite states raise enormous interests in 
quantum information processing and one-way universal 
quantum computing [38]. They are essential for several 
quantum error codes and communication protocols [39], 
as they are robust against decoherence. 

In spite of the previous unbalanced present status be- 
tween entanglement and nonlocality, the relation between 
both quantities for three qubits is seen in a new light 
when we study both resources for those states that at- 
tain the maximum possible nonlocality value given by 
expression (36). One way is to consider what class of par- 
ticular states (35) is maximally nonlocal for a given value 
of their degree of mixture, which is a tool employed to 
characterize mixed states. If we choose the participation 
ratio R = 1/Tr(p 2 ), we can obtain what is the functional 
form of (36) in terms of R. 

This procedure is virtually identical to the variational 
calculation performed for two qubit mixed states. Fol- 
lowing the exact treatment as in Eq. (12), and taking 
into account that Tr[B Mermin ] —32, we obtain 



Mermin max (R) cx 



8 — R 
8R ' 



(38) 



The constant in (38) is obtained by requiring 
Mermin maK to be equal to 4 for pure states (R = 1). One 
class of states that possess the previous optimal value is 
pdiag _ ^ _ 7 x ,x,x,x,x,x,x,x), which is, interestingly 
enough, the generalized Werner state for three qubits 



3=0 



■ n = 3 - x\GHZ)(GHZ\ 



Pw 



(39) 



The exact form for Mermm max is rather unpleasant. 
In practice, the previous bound is an excellent one, 
differing from the exact one by a small amount and 
being equal in those cases where we have a high degree 
of symmetry in the state. For most practical purposes, 
one can consider the equality in (36) to hold. 



where Is is the 8x8 identity matrix, and x = 1 — 8x. 
Notice that this was not the case for the two qubit in- 
stance. 

This interesting feature enables us to discuss the dif- 
ferent ranges for R where to compare nonlocality and 
presence of genuine tripartite entanglement. On the one 
hand, from (38) we obtain the nonlocality critical value 
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i?i = 32/11 ~ 2.9: no three qubit states possess any non- 
locality for participation ratios R > R\. On the other 
hand, the special nature of generalized Werner states al- 
low us to compute the separability threshold between en- 
tanglement and separability [40]. From Ref. [40], the 
contribution x in (39) is such that x < 1/5 involves ab- 
sence of entanglement. Translated into i?-language, it 
implies a second critical value i?2 = 25/4 = 6.25: no 
three qubit states possess entanglement for participation 
ratios R > R\. 

Therefore, the range of i?-values splits into three 
regions: i) between 1 (pure states) and R\, maxi- 
mum amounts of nonlocality imply the presence of 
entanglement; ii) between R\ and i?2 we have no 
violation of the Mermin inequality, yet there exists 
entanglement; finally, iii) region between R 2 and R = 8 
(maximally mixed state) displays absence of both 
magnitudes. Notice, however, that in appearance there 
is some room left for LVM to hold in the second re- 
gion, where no violation of the Mermin inequality occurs. 

The exploration of nonlocality and entanglement for 
three qubit states would be incomplete without an spe- 
cific example of their applicability, such as the possible 
information-theoretic tasks limitations imposed by the 
former on the latter. There is one such case, which is the 
well-known infinite XY model in a transverse magnetic 
field [32]. This instance was explored in detail in Ref. 
[31]. 

The XY model is completely solvable, a fact that al- 
lows us to compute the reduced density matrix for three 
spins without the explicit construction of the global infi- 
nite state of the system. The reduced state of three spins 
reads as 



(a,6) 
Pijk 



(40) 



where i < j < k indicate the positions of the three spins 
and a = j — i,b = k — j their relative distances, {u, v, w} 
denote indexes of the Pauli matrices {tro, o~ x , o~ y , cr z }, and 

Tuv'w = {o~u <8> cr J v (8) a^,}ab- The calculation of the three- 
spin correlations Tuv'w were computed in Ref. [31] by 
using the Wick theorem in quantum field theory. 

The most significant result is that, for any value of 
the anisotropy and external magnetic field, we have 
Mermin™ 13 -* to be less than or equal to 



VVi^) 2 + 4(&) 2 + 4(7&i 6) ) 2 + 4(TiSi 6) ) 2 , 

(41) 

which is always < 2 for any configuration of the spins 
(a,b) yet there is no null entanglement. This constitutes 
a clear sign that states (40) never violate the Mermin in- 
equality, which entails an inherent limitation to the use- 
fulness of entanglement itself. Furthermore, these states 
are shown to be distillable in most of the cases, which is a 



novel result: we have three-party distillable states in the 
XY model with no violation of the Mermin inequality. 
The distillability issue constitutes the subject of the next 
section. 



B. Distillability and nonlocality for three qubit 
states 

Distillability and the violation of Bell inequalities - 
nonlocality- constitute two manifestations of entangle- 
ment. While the former is related to the usefulness in 
quantum information processing tasks, due to the fact 
that most of them require pure-state entanglement as a 
key ingredient, the latter expresses the fact that a state 
cannot be simulated by classical correlations. In this 
vein, Gisin relates both resources when he points out in 
Ref [41] the question of wether there exists any bound 
entangled state that violates some Bell inequality. By 
bound entangled state one implies a state that cannot be 
distilled by means of local operations and classical com- 
munications. In the bipartite case, it has been shown [42] 
that no bound state violates the CHSH inequality. 

The separability criteria borrowed from the bipartite 
case, which employ positive partial transposition [43] for 
all parties, are all approximate. Nevertheless, this cri- 
terion based on the positivity of the ensuing partially 
transposed matrix p Tj has a very interesting application. 
Notice that if a three qubit state p has positive p Tj for 
all j = 1,2,3, where Tj represents partial transposition 
for the system j, then it is said that the state is GHZ- 
distillable, that is, one can distill a GHZ state from many 
copies of p by LOCC [44]. 

One of of questions that we want to address is wether 
there exists any nonlocal bound entangled state of three 
qubits. After applying a series of local transformations, 
one can convert any state into one belonging to the fam- 
ily p^Mermin (35)- For a state of three qubits to be non- 
distillable (bound entangled) , any of the subsequent fol- 
lowing inequalities must hold: 



p Tl > 



^2 + ^2 > A^ — A 1 

A3 + A3 > A,| - A 



P T3 > : 



> 



A 2 + 



>A+ 



A3 + A3 > Af 

'A+ + Ar >A+ 
A3 + A 3 



^0 



> A^ — A 2 



A+ + A+ < 



A+ + A+ < 



A+ + A+ < 



2' 
1 

2' 
1 

2' 
(42) 



where the last inequality in each case is a consequence 
of the sum of the previous two. None of the previous 
inequalities for the eigenvalues of states diagonal in the 
Mermin-basis is compatible with (34) being greater than 
2, which implies that no bound entangled state is present 
in Mermin- diagonal mixed states of three qubits that vi- 
olates the Mermin inequality. 
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1.8 




Mermin max 

FIG. 3. Probability (density) distribution for the bound to 
nonlocality measure Mermin" 1 ^ (34) for all pure and mixed 
states of three qubits diagonal in the Mermin basis. See text 
for details. 



A complementary Monte Carlo numerical survey was 
performed over a set of 10 9 sample states generated 
with a uniform distribution for the set {A^} of the con- 
comitant eigenvalues. This exhaustive, random explo- 
ration confirms the previous result. Fig. 3 depicts the 
probability (density) distribution for nonlocality measure 
Mermin ma * (34). Notice the strong biased behavior to- 
wards no Mermin inequality violation, as well as the rel- 
ative scarcity of those states with some nonlocality. 

A previous work [45] considered too the connexion be- 
tween distillability and violation of the Mermin inequal- 
ity for three qubits. It was concluded there that for a 
particular 4 parameter three qubit states, nonlocality im- 
plied Mermin-distillability. With our analysis, which em- 
braces a more general class of states, we find that this 
is not the case. In point of fact, our numerical explo- 
ration obtains a probability 0.293 to find distillable states 
with no violation of the Mermin inequality, whereas those 
states that being distillable and achieving some nonlocal- 
ity nearly possess a zero-measure (probability of 0.008). 
The vast majority of states are found, with probability 
0.698, both bound entangled and with no violation of the 
Mermin inequality. 

In view of our results, it seems plausible to assume 
that mixed states of three qubits with high amounts 
of Mermin nonlocality, which are likely to be Mermin- 
diagonal, possess no bound entanglement and are thus 
non-distillable. 



IV. FOUR QUBITS 

The first Bell inequality for four qubits was derived by 
Mermin, Ardehali, Bclinskii and Klyshko [46]. It con- 
stitutes of four parties with two dichotomic outcomes 



each, being maximum for the generalized GHZ state 
(|0000) + |llll))/v / 2- The Mermin-Ardehali-Belinskii- 
Klyshko (MABK) inequality reads as Tt(pBmabk) < 4, 
where Bmabk is the MABK operator 

-Sim — -B1112 — -B1121 — -B1211 — -B2111 ^ -B1122 — -B1212 

— -B2112 — #1221 — -B2121 — -B2211 + i?2222 + -B222I 
+ -B2212 + -B2122 + -B1222, 

(43) 

with B uvwx = u-<7®v-cr(g)w-<7®x-<7 with a — (a x ,a y , a z ) 
being the usual Pauli matrices. As in previous instances, 
we shall define the following quantity 



MiBr" = max Tr(pB M AB K ) (44) 

as a measure for the nonlocality content for a given state 
p of four qubits. aj and bj are unit vectors in K 3 . MABK 
inequalities are such that they constitute extensions of 
previous inequalities with the requirement that gener- 
alized GHZ states must maximally violate them. New 
inequalities for four qubits have appeared recently (see 
Rcf. [47]) that possess some other states required for 
optimal violation. In the present study we limit our in- 
terest to the MABK inequality, although new ones could 
be incorporated in order to offer a broader perspective. 
However, with respect to entanglement, little is know for 
the quadripartite case, and thus little comparison can be 
done. 



A. Nonlocality for four qubit states. Extension to 
generalized GHZ states 

The maximization of the MABK inequality 
Tt(pBmabk) < 4 for four qubits mixed states is 
done along similar lines as previously performed for 
bipartite and tripartite cases. Demanding maximum 
amount of nonlocality (44) is tantamount as com- 
puting their optimum values for those states that 
concentrate nonlocal correlations. In the case of four 
qubits, those maximally correlated states are the 
ones which are diagonal in the BelLi basis defined by 
|*±> = ± \l5-j))/V2- Therefore we shall consider 
the following class of four qubits mixed states 

fi K - bawwi + Aji*,- wo, (45) 
3=0 

with ordered eigenvalues Aj + i > Aj. 

Computation of one term B a ^ 7 s = a.-a®f3-<j®^-a®8-a 
of the MABK operator (43) for the BelU basis reads as 
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mixed states is equal or less to 



(46) 



(B aM s)o = a z {3 zlz 5 z ± Re[a+/3+ 1 +S+} 7 
(B a ^s)t = -a z p zlz 5 z ± ite[a + j 8 + 7+5-] ) 
(£a/M>2 = -a z f3 zlz 5 z ± i?e[a + /3+ 7 -5+], 
(B a/ 3 lS )f = a z (3 z j z 5 z ± Re[a + /3 + j~6~], 
(B a p^s)f = -a z /3 z j z S z ± i?e[a _ /3 + 7 _ (5 _ ], 
{B a p~ts)i = a z (3 z j z 5 z ± Re{a+ f3--{ + 5-}, 
(B afh 6)i = a z /3 zlz 5 z ± i?e[ a -/? + 7 + <T], 
(£a/M>7 = -a t MzS, ± ■ReK/3 + 7 + ( 5+], 



Gathering all pure-state expectation values, 
Tr {p { MABK®MABK) is of the form 



7 

^ ±/(n|)(A+ + AT) + <?(a, /3, 7, <J) (At + AJ), (47) 
j=o 

with /(•) is a real function of the product of all z- 
components of the four parties' settings, and g(-) rep- 
resents a real function of several products of all parties' 
components, all of them according to the special form of 
the MABK Bell inequality operator (43). 

The maximum value of (47) is attained with = 0, 

that is, no z-depcndcncy. Similarly to the calculations 
carried out in Appendices I and II, the optimum value of 
the violation of the MABK inequality Tr{pBuABK) < 4 
for mixed states (45) is of the form 



max Tr(pBuABK) ^ 4v^ 



A") 2 - 



(48) 



The exact form for (48), obtained by recourse to convex 
optimization, is extremely complicated (combination of 
rational functions of radicals involving integer powers of 
(At — XT)). As in the case of three qubits, the previous 
bound is an excellent one, and hence we can consider the 
equality in (48) as a close or exact measure of the amount 
of nonlocality present in a BelU diagonal mixed state of 
four qubits. 

When comparing this last result for four qubits 
with those of two and three qubits, we see that all 
three cases involve the same functional form for the 
eigenvalues of the mixed multipartite state diagonal in 
the concomitant maximally correlated basis. This is 
not surprising since the Bell inequalities considered so 
far are multipartite generalizations of the CHSH Bell 
inequality [46]. Therefore we can conjecture the form for 
the maximal violation of the n-party generalized MABK 
inequality Tt(pBmabk) < 4 for a mixed state diagonal 
in the corresponding maximally correlated basis. 

Conjecture. The maximum amount of violation of n- 
party generalized MABK inequalities for diagonal n-qubit 



2^ 



\ E (tf-Aj)' 

\ 3=0 



(49) 



Generalized GHZ states are those states of N parties 



+ y/T^ft"- 1 - 1) = cos a\0) + sin c^" 1 - 1) 

(50) 

who maximally violate the MABK inequality [46] 

stands 



LVM 
MABK,, 



where B L J™ 



Tr{pB M ABK N ) < B 
for the maximum violation allowed by a local variable 
model. We denote the quantal maximum violation [46] 

^ B MABKv =22. 

The maximum value for the corresponding MABK in- 
equality for states (50) is obtained by exactly following 
the optimization procedures carried out in the Appen- 
dices. In point of fact, GHZ states (50) are linear com- 
binations of two pure states that can be written in the 
corresponding maximally correlated basis for that par- 
ticular number of parties. That is, we can rewrite for 
convenience states (50) in the form 




2 V 2 



(51) 



with Ai,2 = \ ± \/p(l — p) (Ai > A2). This new form 
enables us to treat generalized GHZ as linear combina- 
tions of maximally correlated states in a similar fashion 
as performed for two qubit states. 

After some algebra, we obtain that the leading term 
in the violation goes as 2 B m A abk n — P) (sym- 
metric around p = |), from which we reobtain, after 
equating it to B m v a a ^ Kn , the well known result [48] 

sin 2a < l/y/2 N ~ 1 . Thus, by employing our optimiza- 
tion procedure, we not only recover the range where 
generalized GHZ states violate a Bell inequality, but 
also obtain its exact amount. 

As far as entanglement for four qubits is concerned, 
we encounter a considerable discrepancy between max- 
imum entanglement and nonlocality. First of all, it is 
well known that no proper entanglement measure is op- 
erational yet for states living in arbitrary Hilbert spaces. 
However, some measures (for pure states) based on par- 
titions of the system have been advanced, such as the 
so-called global entanglement (GE), which describes the 
average entanglement of each qubit of the system with 
the remaining N — 1 qubits. The GE measure is widely 
regarded as a legitimate A-qubit entanglement measure 
[49-51]. Having a general mixed state implies that the 
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previous measures do not apply as such. To overcome this 
fact we require the extension of these partition-based en- 
tanglement measures by recourse to the usual convex roof 
[52] defined over some given set of pure states. Given the 
extraordinary numerical effort that this procedure would 
imply, an alternative measure is given by the sum of the 
von Neumann entropy of the reduced density matrices of 
all individual qubits, that is, S v n = ^vN(Pi)- Inci- 
dentally, it has also been considered as a proper entan- 
glement measure in the literature [53]. 

It is in this sense that, when employing the sum of 
the von Neumann entropy of all partitions of a state, 
the maximum is reached [54] for a particular (pure) state 
different from the generalized GHZ state for four qubits 
(which is not the case for two or three qubits) . Therefore 
we encounter that maximum entanglement does not cor- 
respond to maximum nonlocality for four qubits already 
at the level of pure states. The analysis for mixed states 
-not performed here- would simply confirm this result. 



V. CONCLUSIONS 

In the present work we have studied how nonlocality is 
present in systems of two, three and four qubits. We have 
exhaustively explored several aspects that are shared by 
quantum entanglement and nonlocality -measured by 
the maximum violation of a Bell inequality- as well as 
pointed out those ones that differentiate these two mag- 
nitudes. By highlighting those issues that concern en- 
tanglement and nonlocality, we shed a new light on the 
connections that exists between these two concepts that 
play a paramount role in quantum-information theory 
and, in turn, in the foundations of quantum mechanics. 

By means of a new optimization method, we have com- 
puted the maximum violation of the CHSH inequality for 
two qubit systems and obtained the concomitant maxi- 
mal states MNMS within the context defined by the de- 
gree of mixture, as measured by either the participation 
ration R of the maximum eigenvalue A max of the state p. 
The direct comparison with MEMS states illustrates an 
anomaly that appears between entanglement and nonlo- 
cality already for mixed states of two qubits, enhanced 
by the information-theoretic tasks limitations that ap- 
pear in the study of bipartite states in the infinite XY 
model. The study of nonlocality for linear combinations 
of pure states of two qubits allowed us to compare how 
both nonlocality and entanglement are distributed in the 
pure state superposition. 

The extension to three qubit states was done along sim- 
ilar lines, employing the maximal violation of the Mermin 
inequality as a nonlocality measure. Analogous compu- 
tations allowed us to obtain the expression for the nonlo- 
cality present in mixed states diagonal in the GHZ basis. 
We obtained that the generalized Werner states for three 
qubits possess maximum nonlocality for a given value of 
the degree of mixture, contrary to the two qubit scenario. 
Also, we extend a previous result concerning distillability 



and nonlocality in the light of quantum entanglement. 

The study of four qubit systems was performed fol- 
lowing the same steps as in the previous two cases: by 
maximizing the MABK inequality for four qubits, we ob- 
tained the maximum violation of this nonlocality measure 
for mixed states diagonal in the BelLi basis. As a conse- 
quence, a careful quantitative analysis is performed for 
generalized GHZ states as well. As far as entanglement 
is concerned, we observe the first discrepancy between 
maximum entangled and optimal nonlocality already for 
pure sates of four qubits. Obviously, the MABK inequal- 
ities are not the only existing Bell inequalities for states 
of arbitrary number of qubits, but because it constitutes 
a simple generalization of the CHSH inequality (never- 
theless, some authors introduce other inequalities that 
incorporate the MABK ones as special cases [55]), it has 
been enough to make use of this particular family of in- 
equalities to illustrate the basic results of the present 
work. Some further work will be required regarding dif- 
ferent multipartite Bell inequalities. 

Despite the fact that for small quantum systems we 
recognize a simple correlation between entanglement and 
nonlocality, the entire situation becomes more involved 
when the dimension of the Hilbert space of the system 
or subsystems augments. This fact is certainly transcen- 
dental for several information-theoretic task require the 
presence of either quantities. Physical situations such as 
the one encountered in the XY model, were null nonlocal- 
ity for two or three parties is compatible with non-zero 
entanglement, do not contribute to unify the ultimate 
quantum correlations that define the state of a quantum 
system. Rather, we are tempted to regard nonlocality 
and entanglement as different quantum resources in view 
of the undefined limits between them. 

On the whole, however, many aspects that also concern 
nonlocality and entanglement have not been considered 
here. One such example could be the so called monogamy 
of entanglement, a fundamental property stating that if 
two quantum systems are maximally correlated (maxi- 
mum entanglement), then they cannot be correlated with 
a third party. This is, for instance, the basis for se- 
cure quantum key distribution based on entanglement 
[7]. The fact that this trade-off also occurs for nonlocal- 
ity [56] in the multipartite case constitutes an issue that 
surely deserves future study. 
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APPENDIX I 

The goal of this appendix is to derive the maximum vi- 
olation of the CHSH inequality (4) for two qubit systems. 
Such endeavor might render somewhat difficult the study 
of the general instance, but this is not the case for there 
is no need to explore the whole space of mixed states 
of two qubits. Since we require Tr(pBcHSii), which is 
a convex function of the two qubit state p, to be maxi- 
mum, it suffices to consider those states that concentrate 
all quantum correlations after the action of a depolariz- 
ing channel [22]. This class of states are, as expected, 
the Bell diagonal states. 

The optimization is taken over the two ob- 
servers' settings {aj,bj}, which are real unit vec- 
tors in M 3 . We choose them to be of the form 
(sm6*fccos<^fc,sin#fcsin(/>fc,cos#fc). With this parame- 
terization, the problem consists in finding the suprc- 
mum of Tr(pBcHSH) over the {k = 1 • • • 8} angles of 
{a x , bi, a 2 , b 2 } that appear in (5). 

The general form entering the Bell operator (5) for one 
single entry is of the kind a ■ a ® /3 • u. Written in the 
computational basis {|0), |1)}, we have 

/ a z /3 z a z /3~ a~f3 z a~/3~ \ 

a z /3 + ~a z f3 z a~/3 + -a~f3 z \ . , 

a+(3 z a+p- ~a z (3 z -a z p~ ' ^ 

\a+/3+ -a+/3 z ~a z f3+ a z /3 z J 

with <x^ = a x ± ia y and ^ = f3 x ± i(3 y . The evaluation 
of (52) for all states in the Bell basis reads as 

(a-<T®0-a)*±=a x p z ±Re[a + + ], 
(a- a® p ■ a)v± = -a z (3 z ± Re[a+ /?"]. 

The expression for Tr(pg l J l 9 ^ Be h s h) , with eigenvalues 
Ai > A 2 > A 3 > Azt, can be cast as 

((Ax - A 4 ) - (A 2 - Aa))[oJ(6? + b x 2 ) + a x 2 (b* - b x 2 )] 
+ ((A 2 + A3) - (Ai + A 4 )) [a\{b\ + b\) + a\{b\ b\)\ 
+ ((Ai - A 4 ) + (A 2 - A 3 )) \a\{b\ + b z 2 ) + a\{b\ - b\)\. 

a z f5 z + Re[a + /3+] ilm[a+ '/?+] 

-Um[a + (3+] a z /3 z - Re[a+ P+] 

-iIm[orp x + a z (3-} Re[a z p+ + a+/3 z } 

Re[a z p+ ~ a+/3 z ] -ilm[a x 0- + a+(3 z ] 



Eigenvalue coefficients in the first and second terms of 
(54) are strictly positive, whereas in the second one it is 
undefined. By rearranging terms in (54), we obtain 



(\ 1 -\ i )[al(bl+bl) + al{b\ -b z 2 ) 
+ (aWi+b x 2 ) + aWi-b x 2 )} 

+ (A 2 - A3) [a\{bl + V 2 ) + a z 2 {b\ - b z 2 ) (55) 
-{al{bl + bl) + aWi-bl))] 

+^[a\{b\ + bl)+al{b\ r bl)], 
with A\ — (A 2 + A3) — (Ai + A4) possessing no clear sign, 
which would imply some further insight. However, on the 
contrary, this fact points out that no y-dependency makes 
(55) even greater. Also, the symmetry in (55) allows us 
to choose the alignment of one of the settings. From 
inspection of (55), we therefore optimize it by choosing 
{ ai = (l,0,0),a 2 = (0,0,-l),b! = (&f,0,-6£),b2 = 
(&f,0,6|)}. 

The final concomitant result amounts to 



max Tr(pBcHSH) = 

max 2(A a - \±)[b z 2 + b x ] + 2(A 2 - A 3 )[&* - b x ] = 
max [2\/2(Ai - A 4 ) cos(6) + 2V2(X 2 - A 3 ) sin(6)] = 
2V2 y /(\ 1 - A 4 ) 2 + (A 2 -A 3 ) 2 . 

(56) 



In the case where nonlocality is to be found in linear 
combinations of pure states of two qubits, we shall per- 
form a similar analysis. Our starting point is the matrix 
of expectation values of elements a ■ a <£> (3 ■ a (52) in the 
Bell basis {|$+), |$"), |*+), |*-)} 



Um[a~l3 z + a z /3~} Re[a z /3 + - a + (3 z ] 

Re[a z j3 + + a + j3 z ] Hm[a z /3~ + a + j3 z ] 
-a z p z + Re[a+(3-] Um[or 

-ilm[a-j3+] -a z (i z - Re[a+j3-\ 




Let us consider a general pure state of the form 

10) = \/A^|$ + >+v^|*~>+V^3l* + >+V^I*~>> (58) 

with real coefficients { -\/Ai} such that Ai + A 2 + A 3 + A 2 = 



1. We will demand the latter to be sorted in decreasing 
value though, as we shall see, this is not mandatory. The 
general case with complex coefficients is somewhat more 
involved. However, for most practical purposes, it will 
suffice to consider real states of the form (58). 
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The fact of having real coefficients in (58) greatly sim- 
plifies the expectation value (<t>\BcnsH\<f>) ■ Its general 
term term {<f>\a ■ a ® (3 ■ a\4>) (52) is of the form 

\i{a z l3 z + Re[a + I3+}) + A 2 (a z /3 Z - Re[a+/3 + ]) 
A 3 ( - a z f3 z + Re{a+f3-}) + A 4 ( - a z p z - Re[a + 
+ ^X~ 3 2Re[a z /3 + + a+(3 z ] + v %X~ 4 2Re[a z l3+ - a+/3 z 

m 



Notice that, in view of (59), that the optimization of the 
CHSH inequality for linear combinations of pure states 
in the Bell basis is almost identical to the corresponding 
mixed state case (differing in the last two terms of (59)). 



Proceeding as before, optimization of ((/)\Bchsh \4>) re- 
turns the value 



2v^/ (Ai - A 4 ) 2 + (A 2 - A 3 ) 2 + 4[v/a7a;] 2 + A[^\tf = 2^2^/^ + A 4 ) 2 + (A 2 + A 3 ) 2 



(60) 



APPENDIX II 

In this appendix we shall derive the explicit form 
for the maximum amount (34) of violation of the Mer- 
min inequality for a three qubit state. As expected, 
since ^(pB Mermin) is a convex function of the quantum 
state p, its maximum is obtained only for pure states, 
namely, the whole class of states forming the Mcrmin- 
basis I**) = ± |7 - j))/V2. In view of this obser- 
vation, we shall consider instead what is the maximum 
violation attained for mixed states diagonal in this basis, 

that is, states of the class p^jermin (35)- Another argu- 
ment for studying these states is that any initial state p 
can be converted into one in the class by means of LOCC. 

Optimization of Mermm max (34) for states PMermin 
(35) is carried out in the same fashion as in the previ- 
ous bipartite case. Once the observers' settings {aj,bj}, 
which are real unit vectors in K 3 , are parameterized in 
spherical coordinates (sn\9k cos (fik, sin 9 k sin cfik, cos 9 k), 
the problem consists in finding the supremum of (34) 
over the set of {k = 1 • • • 12} possible angles for 
{ai,b 1 ,a 2 ,b2,a3,b 3 } in (33). 

To start with, let us write a generic element of the 
Mermin operator (33) of the form B a /3 7 = a-(T<E)/3-a<E)j-a 
in the computational basis {|0),|1)} defined by the z- 
projections of a = (a x ,a y ,a z ). B tt| g 7 reads as 



a z /3 z 
a z (3+ 
a+p z 

+ R+ 



a 



a z (3 

a+/3- 
-a+/3 z 



a (3 Z 
-a z p+ 



with a* = a x ± ia y , (3 d 



crfT \ 

-a j3 z 
-a z p~ 
a z p z J 



± if3 y , and 7 d 
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(61) 



being "rising" and "lowering" terms in the x-y plane. 
Nonlocality measure (34) for diagonal states (35) is com- 
puted by recourse to four B a p 7 s in (61) for different con- 
figurations of vectors. 

The evaluation of (34) for diagonal states (35) amounts 
to compute the expectation value (^flBa/j^l^f) = 

1 _|_ 

(af3^)j for all states {l^ )} in the Mermin basis and 
several vector configurations. The positions in (aftj)^ 1 
are such that a, (3 and 7 correspond to the first, sec- 



ond and third observer, respectively. This computation 
returns 

(a/3f)% = ±i?e[a+/3 + 7 + ], (a^f)f = ±Re[a+P + 1 -] 

(afa)i = ±Re[a+p-j+], (afa)f = ±Re[ a -p + 1 +]. 

(62) 

As we can observe, the observers settings can be two- 
dimensional (he have no z-dependency). 

From previous definitions, we now write 

Tr(p Mermin B Mermin) as 

(A+ - A -)[(aaa)+ - (abb)+ - (bab)+ - (bba)+] 

+(A+ - Ar)[(aaa)+ - (abb)+ - (bab)+ - (bba)+] 

+(A+ - A 2 -)[(aaa)+ - (abb)+ - (bab)+ - (bba)+j 

+(A+ - A3) [(aaa)+ - (abb)+ - (bab)+ - (bba)+] . 

(63) 

The explicit evaluation of the previous quantity can be 
cast as 

( A o + /i - V/i) - «K T «af + a\al)al\ 

-[{albl-a\bl)b%T{atbl + a\bl)bl] 
-[(blal-b\al)b%T{blal + b\al)bl] 

-[(W-M«St(^ + W] 
+ 

( A 2 + /3 - a 2 - /3 )[[«4 + «K ± ««! - 

-[{am + a\b\)bl ± [a\b\ - a\b x 2 )b\\ 

-[{blal + b\al)b%±{btal-b\al)bl] 

-[{blbl + b\bl)a%±{btbl-b\bl)al\]. 

(64) 

Since the Mermin inequality settings are such that it 
must posses rotationally invariance (x-y plane), we are 
free to fix one of them. In view of (64), we choose 
a 2 = (—1,0,0). Also, differences in each term of (64) 
must be maximum in absolute value, which is compat- 
ible with fixing b 2 = (0,1,0). Further calculations im- 
ply a configuration of the type {ai = (of, a\ , 0), a 3 = 
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(af,a|,0),bi = (-of,a?,0),b 3 = (a|,-a|,0)}. Ex- 
pected value Tr{p^^l im B Me rrmn) greatly simplifies 
from (64) into 



( A o + /i - Vi)[2|ofaf | + 2|a?af| ± 2|a?af(l ± 2|« 



(A+3 - A- /3 )[2|<af | - 2|<af | T 2|afa|f| ± 2|o?o|| 



(65) 



Notice that we have reduced our optimization problem 
to one which entails only two real quantities. By intro- 
ducing explicit angles (<j),ip), and after some algebra, we 
obtain 



max Tr(pBMermin) = max (Aq — A )4 sin </> sin 

aj.bj 4>,if> 

+(A+ — A^)4sin</>cos'i/' 
+( A 2 — A 2 ")4cos(/)cosV' 
+ (Ag — A3 )4cos0sin^>. 

(66) 

The solution of (66) is obtained by recourse to the use 
of convex optimization techniques [57] . We do not worry 
about the signs in each term of (64) since we have chosen, 
without loss of generality, the eigenvalues {A*} to be 
sorted in decreasing value. By solving the set of equations 



tan < 



tani/> = 



/ 4(A+-A -)tany> + 4(A+-Ar) 
- Ag ) + 4(A+ - A3 ) tanV> 
/4(A+-A o -)tan0 + 4(A+-A 3 -) 



V4(A+ 



A 2 ) + 4(XJ - X 1 )tan<£ 



(67) 



we finally obtain the desired evaluation of (34) for diag- 
onal states (35). Though the final result is rather cum- 
bersome, we nevertheless derive an excellent bound. In 
view of the coefficients in (66) (the sum of their squared 
values equals one), we provide the final result (36). 
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